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ABSTRACT 

We use data from the Rossi X-ray Timing Explorer to search for harmonics and side- 
bands of the two simultaneous kilohertz quasi-periodic oscillations (kHz QPOs) in Sco 
X-1. We do not detect any of these harmonics or sidebands, with 95% confidence 
upper limits to their power between ~ 1 % and ~ 10 % of the power of the upper kHz 
QPO. The oscillations produced at these frequencies may be attenuated in a scattering 
corona around the neutron star. We find that upper limits to the unattenuated power 
of some of the strongest theoretically predicted harmonics and sidebands are as low 
as ~ 2 % of the unattenuated power of the high-frequency QPO in Sco X-1. 
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1 INTRODUCTION 



It is four years now since the kilohertz quasi-periodic oscil- 
lations (kHz QP Os) were discovered in th e persistent flux 
of Scorpius X-1 Jvan der Klis et al. 1996a|) and 4U 1728-34 
( Strohmayer, Zhang fc Swank 1996j ) . In the meantime, sim- 



ilar kHz QPOs have been seen in some 20 other low-mass 
X-ray binaries (LMXBs; see van der Klis 2000 for a review). 
These QPOs often appear in pairs, with frequencies v\ and 
V2 (yi > Vi) between ~ 400 and ~ 1300 Hz, which in a 
given source can shift by a few hundred Hz, apparently as a 
function of mass accretion rate. 

Most of the models proposed so far assume that one of 
the kHz QPOs reflects the Keplerian orbital motion at some 
preferred radius in the accretion disc (e.g., Miller, Lamb & 
Psaltis 1998; Stella & Vietri 1999; Osherovich & Titarchuk 
1999), but there are other explanations as well (Klein et 
al. 1996a,b; Jernigan, Klein & Arons 2000). In recent dis- 
cussions (Lamb & Miller 1999; Stella 1999; Psaltis 1999) it 
was emphasized that empirical discrimination between two 
of the leading classes of models is possible in principle by 
studying the harmonic and sideband structure of the kHz 
QPOs. Here we concentrate only on these two model classes 

In the 'sonic-point' model (SPM; Miller et al. 1998), the 
QPO at v-2 (the upper QPO) is produced at the radius where 
the radial flow velocity in the disc turns from subsonic to 
supersonic (the sonic radius), and the QPO at v\ (the lower 
QPO) originates by a beat between the upper QPO and 
the spin frequency of the neutron star. In the 'relativistic- 
precession' model (RPM; Stella & Vietri 1999) the QPO at 



v% is also assumed to be Keplerian, but the QPO at v\ is 
produced by the apsidal precession of a slightly non-circular 
inner accretion disc. The QPO frequencies in the RPM are 
calculated for test particles in purely geodesic relativistic 
motion, i.e., neglecting the hydrodynamical and radiative 
effects of the accretion flow. However, Psaltis & Norman 
( |2000| ) have recently proposed a dynamical model in which 
the QPOs are produced by oscillations in the accretion disk. 
In this model, which we will call 'transition-radius' model 
(TRM) , there is a transition radius in the accretion disc that 
acts as a band-pass filter with resonances near the orbital 
and periastron-precession frequencies. 

Besides the main peaks at v\ and v%, the SPM and the 
TRM predict other (weaker) harmonics and sidebands of 
these QPOs, at specific frequencies. For instance, the SPM 
predicts a relatively strong harmonic of the lower QPO at 
2v\ (see Table 3 of Miller et al. 1998 for a list of other side- 
band peaks predicted by the SPM), whereas the TRM pre- 
dicts a sideband at 2^2 — v\ (cf. eq. [29] in Psaltis & Norman 
2000). In principle, the detection of a QPO at 2^i and a non- 
detection of a QPO at 2v-z — V\ would tend to rule out the 
TRM, whereas the detection of a QPO at 2vi — v\ and a 
non-d etection of a QPO at 2v\ would tend to rule out the 
SPM ( [Miller 2000| ). 

In this Letter, we use data from the Rossi X-ray Timing 
Explorer (RXTE) to search for the predicted harmonics and 
sidebands in the power spectrum of the kHz QPO source Sco 
X-1. Because of its very high flux, Sco X-1 has extremely 
significant kHz QPOs, and hence a very sensitive study of 
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any harmonic structure is possible. We do not detect any 
of these secondary QPOs. The fact that we detect the kHz 
QPOs at v\ and V2 , but none of these other peaks sets severe 
constraints on the models currently proposed to explain the 
kHz QPOs in LMXBs. 



2 OBSERVATIONS 

We used data from the Proportional Counter Array (PCA; 
Jahoda et al. 1996) on board RXTE (Bradt, Rothschild & 
Swank 1993) taken on 1996 February 14, 18, 19, 1996 May 
24 - 28, 1997 March 15, 1997 April 18 - 24, 1997 August 22, 
1998 January 2 - 8, 1998 February 27, 28, 1998 May 30, 31, 
1998 June 1, 2, 1998 July 2 - 5, 1999 January 6, 8 - 11 and 
13 - 16. Each observation (i.e., each part of the data with a 
unique RXTE ID number) consists of data blocks of ~ 60 s 
to ~ 3, 700 s interrupted by passages of the satellite through 
the South Atlantic Anomaly and occultations of the source 
by the Earth. The total usable time was ~ 630 ks. 

To avoid detector safety triggers, telemetry saturation, 
and to reduce the dead-time effects produced by the high 
count rate of Sco X-l, some observations were carried out 
with the source slightly off-axis, with some of the five pro- 
portional counter units of the PCA switched off, recording 
only photons detected by the upper anode chain of the PCA, 
recording only photons from a limited energy range, or us- 
ing a combination of these constraints. In all cases, high-time 
resolution data were available with a time resolution of 0.25 
ms or better. For our analysis below we combined all the 
available data, irrespective of whether they were collected 
using any of the above observational constraints. When sin- 
gle and double-event data were recorded in parallel (see van 
der Klis et al. 1996b), we combined them off-line to enhance 
the sensitivity. 



3 ANALYSIS AND RESULTS 

We divided the high time resolution data into 16 s segments, 
and produced a power spectrum for each of these segments 
up to a Nyquist frequency of 2048 Hz (for ~ 50 % of the data 
we also produced power spectra up to a Nyquist frequency 
of 4096 Hz). For 69% of the observations the power spectra 
were calculated using the full PCA energy band. For the 
rest of the power spectra we used data from selected energy 
bands (24 % of the power spectra were calculated between 
5-18 keV, 5 % between 5-60 keV, and 2 % between 2 - 
18 keV) because we found that the kHz QPOs were more 
significant in those energy bands, or because those were the 
only energy bands available. Finally, we averaged together 
groups of 8 contiguous 16 s power spectra to produce average 
power spectra. 

We searched these average power spectra for kHz QPOs, 
at frequencies >250 Hz. We did this by first identifying those 
power spectra that showed a strong QPO, and visually es- 
timating its frequency, 1/0 ■ When two QPOs were present in 
the power spectra, we always picked the one at higher fre- 
quency. It turned out that in the cases where only one kHz 
QPO was visible, it was always the upper kHz QPO (see be- 
low). We then fitted the power spectra, in the range vo — 100 
Hz to vo + 100 Hz, using a function consisting of a constant, 
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Figure 1. The frequency separation between the kHz QPOs in 
Sco X-l as a function of the upper kHz QPO frequency. 



a power law and one Lorentzian. We discarded the power 
spectra for which the QPOs were less than 3<r significant. 
We note that in this manner we might have discarded data 
with weak QPOs that could have been detected averaging 
together more data. There were 1,384 average power spectra 
with significant kHz QPOs, equivalent to 177,152 s of data. 

For those observations where we detected only one QPO 
in the aver age power spectra we applied the 'shift-and-add' 
technique (VIendez et al. 1998a) to try and detect the sec- 
ond QPO: Based on the frequency of the detected QPO, we 
aligned and averaged all the power spectra of a single obser- 
vation. In those cases, this procedure revealed a second kHz 
QPO, which in all cases was at a lower frequency than the 
QPO peak that we used to align the power spectra, show- 
ing that all the initial frequency measurements in the 1,384 
average power spectra corresponded to the upper kHz QPO. 

Next, we measured the frequency separation between 
the kHz QPOs as a function of the upper QPO frequency, 
V2- We proceeded as follows: We aligned the 1,384 power 
spectra using the upper kHz QPO as a reference; we grouped 
the data in 49 sets, each of them containing ~ 40 to ~ 210 
power spectra, such that V2 did not vary by more than 5—10 
Hz within each set, and we combined these aligned spectra to 
produce an average power spectrum for each set. We fitted 
these 49 power spectra in the range 400 — 1300 Hz using 
a function consisting of a constant, a power law and two 
Lorentzians. The fits were good, with reduced \ 2 < !•! f° r 
279 degrees of freedom, and the significance of both peaks 
was always > 3o\ 

Fig. [I] shows that Au = V2 — v\ , the frequency difference 
between the two QPOs, decreases from ~ 310 Hz to ~ 240 
Hz, as h*2 increases from ~ 840 Hz to ~ 1100 Hz. This figure 
can be compared to fig. 3a of van der Klis et al. (1997), who 



first reported the decrease in Au with V2 in Sco X-l. Because 
we included more data, and because we used the shift-and- 
add technique to measure Av, the errors are smaller in the 
present figure than in the previous one, and some structure, 
particularly between 900 and 1020 Hz, becomes apparent 
(perhaps this structure is related to the 'bump' seen in the 
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Figure 2. Shifted-and-averaged power spectra of Sco X-l. The 
average power spectrum shown in the left panel was computed 
by first aligning the individual power spectra such that a QPO 
peak at 2v\ would always end up at 1320 Hz, indicated by the 
arrow. In the right panel a QPO peak at 2v2 — Vx would similarly 
end up at 1284 Hz (arrow). In both panels, the vertical lines show 
the frequency range used to fit the data (see text). Notice that 
the frequency alignment also altered the apparent shapes and 
frequencies of the strong QPO peaks. 



plot of Av vs. vx of 4U 1608-52 at vx ~ 700 Hz; see fig. 3 in 
Mendez et al. 1998b). 

From Fig. |l| we can read off v\ as a function of U2 ; be- 
cause we already know v-z for each average power spectrum, 
we can calculate the expected frequencies of hypothetical 
signals at multiples of v\ and V2, or at frequencies that are 
combinations of these two frequencies, for each power spec- 
trum. As we described in §[l], some of these frequencies are 
important in the context of the models proposed to explain 
the kHz QPOs. 

We calculated the following quantities: 2vx, 2i>2, V2 — V1, 
u 2 + ui, 2v\ — U2, 2v2 — vx, 2vx + V2, 2v2 + vx and 2(v2 — vx). 
We then aligned, in turn, the power spectra on each of the 
above frequencies, and averaged them together to try and 
detect a signal present at any of these frequencies. We fit- 
ted these shifted and averaged power spectra in a segment 
of ~ 400 Hz that included the expected frequency, using a 
constant plus a Lorentzian with the centroid fixed in each 
case at the expected frequency, and with a fixed FWHM of 
200 Hz, approximat ely equal to the sum of th e FWHM of 
each of the QPOs (van der Klis et al. 1996b). This value 



is about two times larger than the largest FWHM mea- 
sured in the QPOs of this and other sources, and yields 
conservative upper limits to the unobserved harmonics and 
sidebands; tighter upper limits could be obtained by fixing 
the FWHM to a smaller value. No detections result from 
choosing a smaller FWHM in the current data set. In a few 
cases we added a power law to the fit to take account of 
a sloping continuum. The shape of the power spect rum at 
high fr eque ncies is dominated by dead-time effects (Zhang 
et al. 1991:) which, in the case of Sco X-l at a count rate 
that exceeds ~ 25,000 c s" 1 PCU" 1 , are large and not 
yet sufficiently well understood to predict the shape of the 
high-frequency part of the power spectrum accurately. We 
note that a better knowledge of this effect could yield tighter 
upper limits than those that we report in this paper. 

We did not detect a significant QPO at any of the above 
frequencies. As an example, in Fig. ^| we show two power 



Table 1. Upper limits to the observed amplitudes of the harmon- 
ics and sidebands of the kHz QPOs in Sco X-l 

Frequency Predicted frequency range Upper limit 



2vx 


1090 - 


1690 Hz 


0.12 


2V2 


1690 - 


2170 Hz 


0.13 


V2 - Vx 


240 - 


300 Hz 


0.30 


V2 + vx 


1390 - 


1930 Hz 


0.08 


2vx - v 2 


240 - 


600 Hz 


0.26 


2V2 — vx 


1140 - 


1330 Hz 


0.12 


2vx + v 2 


1930 - 


2770 Hz 


0.10 


2v2 + vx 


2230 - 


3020 Hz 


0.10 


2{v 2 -vx) 


480 - 


600 Hz 


0.21 



a 95 % confidence upper limits to the rms fractional amplitude 
of the QPO at the indicated frequency, in units of the observed 
rms fractional amplitude of the QP Q at vo. This QPO has- an 



amplitude between 0.6 % and 2.5 % ( yan dcr Klis ct al. 1997| ) 



spectra that were shifted to 2vx and to 2v2 — vx, the fre- 
quencies of the main high-frequency peaks (besid es the kHz 
QPOs) predicted b y the models of Miller et al. ( |l99g| ) and 
Psaltis & Norman (2000), respectively. 



Upper limits are shown in Table [j]. To calculate these 
upper limits we varied the amplitude of the Lorentzian un- 
til the x 2 °f the fit increased by 2.71 with respect to the 
best-fitting value (95 % confidence level for a single param- 
eter); the upper limits we quote in Table ^ represent the 
change in the amplitude for this change in the \ 2 ■ Because 
of the observational constraints described in and because 
we combined power spectra computed from data in different 
energy bands, we cann ot give upper limit s in units of frac- 
tional rms amplitude (]yan der Klis 1995 ); instead, all the 
upper limits in Table |l| are in units of the observed frac- 
tional rms amplitude of the QPO at V2- Any peaks at the 
predicted frequencies are between 11 and 156 times weaker 
than the peak at V2, which translates into rms amplitude 
ratios between 0.08 and 0.30. 



4 DISCUSSION 

We have measured the frequencies of the two simultaneous 
kHz QPOs, Vx and V2, in Sco X-l in a large data set. We have 
used these measurements to calculate the expected frequen- 
cies of hypothetical signals at harmonics and sidebands of 
these QPOs, as predicted by two classes of kHz QPO models. 
We used a sensitive technique to search for these harmonics 
and sidebands in the power spectra of Sco X-l, but none 
was detected. Their power is between 10 and 100 times less 
than the power of the upper kHz QPO (95 % confidence). 

Oscillations produced close to the surface of the neutron 
star can be attenuated in a scattering corona. The attenu- 
ation is generally larger for oscillations produced by a pen- 
cil beam sweeping the surroundings (beaming oscillations), 
than for oscillations due to actual changes of the luminosity 
with time (luminosity oscillations; Brainerd & Lamb 1987; 
Kylafis & Klimis 1987; Kylafis & Phinney 1989; Miller et 
al. 1988). The attenuation factor (the ratio of the amplitude 
A x of the oscillations at infinity to their original amplitude 
Ao) as a function of frequency v, for a scattering corona of 
a radius R and optical depth r is given by: 

AacMm 2 3/2 xe- x + e~ T (1) 



-Ao.lun 
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[H] 

Table 2. Upper limits to the unattenuated amplitude of the harmonics and sidebands of the kHz QPOs in Sco X-l 



J. 1 h\_.iil. y 


T a 




T Tnrifir limii - 




R fkml" 


Vj LJLJt/l 111I11L 




< 1 


336 


0.15 


8.7 


129 


0.30 


2v 2 


(5.5 


202 


0.27 


9.8 


109 


0.93 


V2 - V\ 


< 1 


< 10 


0.30 


< 1 


< 10 


0.30 


V2 + VI 


5.7 


262 


0.11 


9.1 


121 


0.35 


2v\ - V2 


< 1 


< 10 


0.26 


< 1 


< 10 


0.26 


2u 2 - v\ 


< 1 


352 


0.14 


8.3 


138 


0.23 


2ui + u 2 


7.3 


167 


0.37 


9.8 


109 


1.43 


2u 2 + ui 


7.8 


151 


0.54 


10.2 


103 


2.16 


2{V2 - v\) 


< 1 


< 10 


0.21 


< 1 


< 10 


0.21 



a Optical depth and radius of a possible scattering cloud around Sco X-l, obtained by minimizing eq. |3j at each frequency (see text). 

These are the values that were used to calculate the corresponding upper limits shown in this Table. 
b Upper limits for the unattenuated amplitude of the QPO at the indicated frequency for a luminosity oscillation (see eq. Jlj ) , in units 

of the unattenuated amplitude of the QPO at v%. 
c Upper limits for the unattenuated amplitude of the QPO at the indicated frequency for a beaming oscillation (see eq. |^]), in units of 

the unattenuated amplitude of the QPO at U2- 
In both cases the QPO at U2 is assumed to be a beaming oscillation (Miller et al. 1998) 



for a luminosity oscillation, and by 



Ao 



2 5/2 ; 



An 



1+T 



+ e 



(2) 



for a beaming oscillation, with x = (SnuRr / 'c) 1//2 (Kylafit 



Phiniey 1989) 



We can use these relations to calculate the largest 
unattenuated amplitudes allowed by the data in Table Q 
Because at each frequency the upper limit on the ratio 
of observed amplitudes is Aoo.yt^Moo.bcam (i/ 2 ) (Table |), 
the upper limit on the ratio of unattenuated amplitudes 
^o,y(^)/^4o,beam(^2), can be obtained by minimizing 



Ax,,y(i/)/.40,y(i/) 

,bcam O 2 )/A ,beamO2) 



(3) 



with respect to r and R, with the additional constraint that 
Aoo, beam (V2)/A), beam (V2) > 0.01 (because the QPO at u 2 is 
detected at an amplitude >1 %; van der Klis et al. 1997). 
In eq. H, y stands either for luminosity or beaming (see eq. 
jlj and ||), and we have assumed that the QPO at V2 is a 
beaming oscillation (Miller et al. 1998). 

The upper limits to the unattenuated amplitudes of the 
sidebands and harmonics of the kHz QPOs, both for lumi- 
nosity and beaming oscillations, are shown in Table |^. In 
units of the amplitude of the QPO at V2, the amplitudes of 
these harmonics and sidebands range from 0.11 to 0.54 for a 
luminosity oscillation, and from 0.21 to 2.16 for a beaming 
oscillation (assuming beaming for the QPO at 1/2 )• From the 
minimization of a y {v) at the frequencies listed in Table |^, we 
get that optical depths r ~ 1 — 10 and radii R ~ 1 — 350 km 
for the putative scattering cloud in Sco X-l would be most 
effective at suppressing the sidebands and harmonics. These 
val ues are com parable to those obtained by by Vaughan et 
al. ( |1997| |1998[ ), r < 5 and R < 80 km, for 4U 1608-52 from 
measurements of the time lags of the kHz QPOs (see also 
Kaaret et al. 1999) . We want to mention that the upper lim- 
its given in Table have been calculated assuming properties 
of the scattering corona that maximize the attenuation at 
each frequency independently. In reality, if the corona exists 
it would have fixed values of r and R, so that if it is effective 
at suppressing the oscillations at some frequencies, it will be 
less effective at others. 



In conclusion, there is a remarkable lack of harmonic 
and sideband structure in the kHz QPOs of Sco X-l. None of 
the frequencies predicted by the two main classes of models 
are detected, with upper limits implying they are one to two 
orders of magnitude weaker in power than the main peaks. 
This sets severe constraints on either the purity of the kHz 
QPO generation mechanism in Sco X-l, or the parameters 
of the scattering corona that surrounds the neutron star. 
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